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THE  EFFECTS  OF  ECCENTRICITIES  ON  THE  FRACTURE 
OF  OFF-AXIS  FIBER  COMPOSITES 
by  C.  C.  Chamis*  and  J.  H.  Sinclair* 

National  Aeronautics  and  Space  Administration 
Lewis  Research  Center 
Cleveland,  Ohio  44135 

ABSTRACT 

Finite  element  analyses  were  performed  to  investigate  theo¬ 
retically  the  effects  of  in-plane  and  out-of-plane  eccentricities, 
bending  or  twisting,  and  thickness  nonuniformity  on  the  axial  stress 
and  strain  variations  across  the  width  of  off-axis  specimens.  The 
results  are  compared  with  measured  data  and  are  also  used  to 
assess  the  effects  of  these  eccentricities  on  the  fracture  stress  of 
off-axis  fiber  composites.  Guidelines  for  detecting  and  minimiz- 
ii^  the  presence  of  eccentricities  are  described. 

INTRODUCTION 

Off-axis  tensile  data  for  unidirectional  composites  are  of  con¬ 
siderable  interest  to  the  fiber  composite  community  for  several  im¬ 
portant  reasons.  Some  of  these  are:  (1)  determination  of  the  varia¬ 
tion  of  elastic  properties  and  fracture  stress  (strain)  as  a  function 
of  load  angle  (angle  between  fiber  and  load  directions),  (2)  verifica¬ 
tion  of  composite  macromechanics  theories  for  elastic  properties 
and  for  combined- stress  fracture,  and  (3)  generation  of  fundamental 
information  for  assessing  angleplied  laminate  mechanical  behavior. 


♦Aerospace  and  Materials  Engineer,  respectively,  Composites 
and  Structures  Branch,  NASA. 
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An  extensive  investigation  was  conducted  to  study  the  mechanical 
behavior  and  fracture  characteristics  of  off-axis  fiber  composites 
(refs.  1  and  2).  Specimens  of  MODI/E  (graphite  reinforced  epoxy) 
were  proposed  and  tested  with  load  angles  varying  from  0  to  90°  with 
the  fiber  direction.  It  was  found  in  that  investigation  that  in-plane  and 
out -of -plane  eccentricities  affect  the  stress  (strain)  variation  across 
the  width  of  the  specimen  significantly.  The  effects  of  these  eccentri¬ 
cities  on  the  tensile  strain  (stress)  variation  and  fracture  of  off-axis 
fiber  composites  are  examined  herein  using  finite  element  analyses. 
Also  the  effects  of  nonuniform  specimen  thickness  on  strain  (stress) 
variation^  and  fractures  are  examined  using  finite  element  analysis. 
The  results  obtained  are  compared  with  measured  data  and  are 
used  to:  (1)  assess  the  effects  of  eccentricities  on  off-axis  composite 
tensile  fracture  and  (2)  establish  guidelines  for  detecting  and  mini¬ 
mizing  eccentricities  during  testing. 

m-PLANE  BENDING  EFFECTS 


Off-axis  tensile  specimens  will  tend  to  undergo  in-plane  bending. 
This  is  caused  by  the  coupling  between  normal  and  shear  deforma¬ 
tions:  this  coupling  will  tend  to  deform  the  specimen  in  shear.  How¬ 


ever,  the  grips  prevent  the  specimen  ends  from  shearing,  thereby  in - 

ducing  in-plane  bending.  This  in-plane  bending  induces  axial  stress 

□ 

and  strain  variations  across  the  specimen  width.  These  variations  .  n. 

are  determined  theoretically  herein  using  finite -element  analysis. 

The  finite  element  used  in  the  analysis  is  a  second-order  triangu -  ^  '  iSjj’" 

— — j 

lar  plate  finite  element  with  six  nodes  and  two  displacement  degrees 
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of  freedom  per  node.  A  schematic  of  the  finite- element  representation 
is  shown  in  figure  1.  The  dimensions  used  in  the  analysis  were  those 
of  the  actual  test  specimens.  Those  shown  in  the  schematic  are  for 
the  10°  off-axis  test  specimen.  Note  that  the  finite -element  repre¬ 
sentation  includes  the  tapered  end-tab  portions  projecting  beyond  the 
grip  ends.  Note  also  that  the  finite- element  representation  consists 
of  288  elements,  657  nodes,  and  1314  degrees  of  freedom. 

Finite-element  analysis  results  for  the  axial  stress  variation, 
near  the  end  tab  (node  line  73  to  81,  fig.  1)  are  summarized  graph¬ 
ically  in  figure  2.  These  stress  variations  were  determined  using 
the  fracture  loads  of  the  specimens  and  the  elastic  constants  sum¬ 
marized  in  table  I.  As  can  be  seen  in  figure  2,  the  most  significant 

axial  stress  variation  is  for  the  10°  off-axis  specimen  with  a  maxi- 

Q  9 

mum  difference  of  16.6x10  N/cm  (24  ksi)  from  edge-to-edge 
(46x10^  to  30x10^  N/cm^  (67  to  43  ksi)).  Additional  discussion  on 
this  variation  is  given  in  reference  3.  The  next  most  significant 
axial  stress  variation  is  that  for  the  15°  off-axis  specimen  with  a 
maximum  difference  of  7.8x10^  N/cm^  (13  ksi)  from  edge-to-edge 
(26X10^  to  17X10^  N/cm^  (38  to  25  ksi)).  The  axial  stress  variation 
for  the  remaining  specimens  is  relatively  mild  and  may  be  considered 
as  insignificant.  An  interesting  result  shown  in  figure  2  is  the  stress 
reversal  trend  from  5°  (increasing  left  to  right)  to  10°  (decreasing). 

The  important  observation  from  the  preceding  discussion  is  that 
off-axis  tensile  specimens  show  high  axial  stresses  at  the  edges  near 
the  grips  in  the  5°  to  15°  load  angle  range  and,  consequently,  fracture 
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should  initiate  in  this  region. 

Corresponding  results  for  the  axial  stress  variation  are  shown 
in  figure  3.  Here,  again,  the  significant  axial  strain  variation  across 
the  specimen  width  is  for  the  10°  and  15°  off-axis  specimens.  These 
results  illustrate  the  importance  of  placing  strain  gages  as  close  to 
the  e(^e  as  possible  near  the  end-tab  region. 

Finite-element  results  for  the  axial  stress  variation  at  the 

Specimen  midlength  (center)  are  shown  in  figure  4.  Only  the  10°  off- 

3  2 

axis  specimen  shows  a  significant  variation  (about  9x10  N/cm 
(13  ksi))  from  edge-to-edge  (39x10^  N/cm^  to  30x10^  N/cm^  (56  to 
43  ksi)).  Corresponding  results  for  axial  strain  are  shown  in  fig¬ 
ure  5.  As  can  be  seen  in  this  figure,  only  the  10°  and  15°  off-axis 
specimens  show  significant  variations  from  edge-to-edge. 

The  important  observation  here  is  that  in-plane  bendiig  produces 
significant  axial  stress  variation  at  midlength  only  in  the  10°  off- 
axis  specimen.  The  significance  of  this  observation  is  that  the  P/A 
(fracture  load/cross-section-area)  stress  is  a  very  good  approxima¬ 
tion  to  the  actual  axial  stress  at  the  center  of  the  off-axis  specimens. 
And,  in  addition,  the  fracture  stress  determined  from  P/A  would 
probably  be  on  the  conservative  side.  It  is  important  to  keep  in  mind 
that  these  comments  apply  to  specimens  with  the  gage  length-to- 
width  ratios  tested  herein,  which  were  14  or  greater. 

Comparison  of  finite -element  predicted  axial  strains  with  meas¬ 
ured  data  near  the  specimen  end  tab  at  fracture  load  are  shown  in 
figure  6.  Correspondiig  results  at  the  specimen  midlength  are  shown 
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in  figure  7.  As  can  be  seen  from  these  figures,  the  agreement  is 
reasonably  good  for  the  three  specimens  near  the  end  tab  and  the  60° 
specimen  at  midlength.  However,  the  agreement  for  the  10°  and  30° 
specimens  at  midlength  is  relatively  poor.  The  predicted  results 
are  about  10  to  20  percent  higher  than  the  measured  data  at  the  left 
edge  and  center  and  are  less  than  10  percent  at  the  right  edge. 

Some  factors  that  may  have  contributed  to  this  poor  agreement 
between  predicted  and  measured  fracture  strains  at  midlength  of 
the  10°  and  30°  off-axis  specimens  are: 

(1)  Inability  to  simulate  mathematically  exactly  the  physical 

boundary  conditions 

(2)  Nonlinear  material  behavior  near  fracture 

(3)  Out-of-plane  eccentricities  -  bending  and/or  twisting 

(4)  Variation  in  specimen  thickness 

Item  (1)  was  extensively  studied  via  sensitivity  analysis  in  ref¬ 
erence  3  and  found  to  have  an  effect  of  less  than  5  percent.  Item  (2) 
is  not  believed  to  have  any  significant  contribution  because  the  stress 
strain  curves  (figs.  6  and  8,  from  ref.  1,  Part  I)  are  linear  to  frac¬ 
ture.  Items  (3)  and  (4)  were  investigated  herein  and  are  described 
in  the  next  section.  Note  that  item  (3)  was  also  discussed  in  refer¬ 
ence  4. 

OUT-OF-PLANE  BENDING  AND  TWISTING  EFFECTS 

The  effects  of  out-of -plane  bending  and  twisting  on  axial  strain 
were  evaluated  for  the  10°  and  30°  off-axis  specimens  using  NASTRAN 
(NASA  Structural  Analysis  Finite  Element  Computer  program, 
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ref.  5).  The  NASTRAN  model  of  the  specimen  is  shown  in  figure  8. 
The  NASTRAN  model  consisted  of  657  nodes  (1971  degrees  of  free¬ 
dom)  and  576  quadrilateral  plate  bending  elements,  which  included 
the  tapered  portion  of  the  reinforcing  end  tabs.  Note  that  the  finite- 
element  representation  includes  two  groups  of  elements.  At  each  end 
the  elements  are  0. 159  centimeter  (0. 0625  in. )  long;  these  represent 
the  tapered  portion  of  the  reinforcing  tabs  and  the  first  quarter  inch 
segment  of  the  test  section,  which  is  the  site  of  the  top  strain  gages. 
The  remaining  elements  of  the  representation  are  0.  318  centimeters 
(0. 125  in. )  long.  All  elements  for  this  model  are  0. 159  centimeters 
(0. 0625  in. )  wide.  The  element  size  was  made  small  enough  to  study 
the  zones  where  the  strain  gages  were  located  on  the  actual  specimen. 
The  material  properties  required  for  NASTRAN  were  generated  from 
the  elastic  constants  in  table  1.  The  load  for  both  the  out-of-plane 
bending  and  twisting  moments  was  11.  3  newton-meters  (100  in.  -lb). 
The  value  of  11.  3  newton-meters  (100  in.  -lb)  was  selected  mainly  for 
convenience.  It  corresponds  roughly  to  an  eccentricity  of  a  laminate 
thickness.  The  effects  of  smaller  eccentricities  are  readily  obtained 
by  direct  proportion  since  a  linear  stress  analysis  was  performed. 

NASTRAN  undeformed  and  deformed  plots  due  to  out-of -plane 

o 

bending  moments  are  shown  in  figure  9  for  the  10  off-axis  specimen 
and  in  figure  10  for  the  30®  off-axis  specimen.  As  can  be  seen  in 
these  plots  the  deformation  for  both  bending  and  twisting  are  con¬ 
siderable. 

The  axial  strain  variation  due  to  bending  moments  across  the 
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specimen  width  predicted  using  NASTRAN  is  shown  in  figure  11  (solid 
lines  for  the  10°  and  interrupted  lines  for  the  30°  off-axis  specimens). 
Corresponding  results  for  axial  strain  variation  at  midlength  are 
shown  in  figure  12.  The  curves  in  these  figures  show  that  the  axial 
strain  variation  can  be  significant  near  the  grips  for  both  bending  and 
twisting  and  at  midlength  for  bending.  This  would  tend  to  explain  the 
differences  between  predicted  and  measured  data  shown  in  figures  6 
and  7  and  discussed  previously. 

Thus  we  see  that  out-of -plane  eccentricities  can  contribute  sig¬ 
nificantly  to  the  axial  strains.  Therefore,  care  should  be  taken  to  keep 
them  to  an  absolute  minimum  during  testing  of  off-axis  specimens. 

THICKNESS  VARIATION  EFFECTS 

The  effects  of  specimen  thickness  variation  on  the  axial  strain 
were  investigated  usii^  NASTRAN  and  actual  measured  thickness 
variations  of  the  specimen  (0.  15  to  0. 14  cm  (0. 059  to  0. 055  in. )). 

The  finite -element  model  used  is  shown  in  figure  8  and  has  already 
been  described.  The  results  obtained  for  the  5°  off-axis  specimen 
are  compared  with  those  for  uniform  thickness  in  figure  13.  As  can 
be  observed  from  the  curves  in  this  figure  the  thickness  variation  ef¬ 
fects  are  negligible. 

GUIDELINES  FOR  DETECTING  AND 
MINIMIZING  ECCENTRICITIES 

The  following  guidelines  may  be  helpful  in  instrumenting  speci¬ 
mens  to  detect  the  presence  of  out-of-plane  eccentricities  durii^ 
testing: 
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(1)  For  out-of -plane  bending,  place  strain  gages  back-to-back  at 
the  specimen  edge  (fig.  14). 

(2)  For  out-of -plane  twisting,  place  strain  gages  at  both  edges  on 
the  same  surface  of  the  specimen  near  the  end  tab  (fig.  14). 

(3)  If,  durii^  testing,  the  differences  in  the  readings  from  the 
pair  of  strain  gages  in  (1)  or  (2)  or  both  become  excessively  high 
(say,  more  than  15  percent),  then  stop  the  test  and  realine  the  speci¬ 
men  to  minimize  the  out-of -plane  eccentricities.  The  strains  already 
recorded  can  be  used  to  guide  the  direction  of  the  realinement. 

SUMMARY  OF  RESULTS  AND  CONCLUSIONS 

The  major  results  and  conclusions  of  this  investigation  are; 

1.  The  second-order  triangular  finite -element  predicted  results 
showed  that  in-plane  bending  has  considerable  influence  bn  the  axial 
strain  variation  across  the  width  of  the  specimen.  This  influence  is 
most  significant  in  the  5*^  to  30®  load-angle  range.  The  predicted 
fracture  strain  variation  was  off  by  about  20  percent  from  the  meas¬ 
ured  data. 

2.  NASTRAN  predicted  results  showed  that  thickness  variations 
in  the  specimen  (0. 14  to  0. 15  cm  (0. 055  to  0.  059  in. ))  have  negligible 
effect  on  the  axial  strain  variation  across  the  specimen  width. 

3.  NASTRAN  predicted  results  showed  that  out-of-plane  bending 
and  twisting  eccentricities  have  significant  effects  on  the  axial  strain 
variation  across  the  width  for  specimens  in  the  10®  to  30®  load-angle 
range. 

4.  Care  should  be  taken  to  minimize  eccentricities  that  will  in¬ 
duce  out-of-plane  bending  and  twisting  since  these  eccentricities 
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have  significant  effect  on  the  axial  strain. 

5.  Fracture  stress  of  off-axis  tensile  specimens  determined  by 
load  to  area  ratio  should  be  on  the  conservative  side. 
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TABLE  I.  -  PREDICTED  COMPOSITE  ELASTIC  CONSTANTS  -  STRUCTURAL  AXES  FOR  MOD  I/E 

[Used  in  finite -element  analyses.  ] 

■; - - ' — I - - - - - — - - 

specimen  Load  Composite  elastic  constants 

angle,  - - - - - 

deg  _ _ _ Poisson’s  Couplii 

p,  ratio,  coefficie 

- p _ _ ^cxy  - 1- 

_ N/cm^  psi  N/cm^  psi  N/cm^  psi  ""cxs  ^ 

0  0  21.0X10®  30.4X10®  0.73X10®  1.06X10®  0.515x10®  0.747X10®  0. 260  1 0  In 

5  16.3  23.6  .73  1.06  .519  .753 

10  9.79  14.2  .74  1.08  .531  .770 

15  5.97  8.66  .77  1.11  .550  .798 

30  2.12  3.07  .88  1.28  .638  .926 

45  1.20  1.74  1.20  1.74  .696  1.01 

60  .88  1.28  2.12  3.07  qoa 

.77 
.73 


30.4X10' 

23.6 

14.2 

8.66 

3.07 

1.74 

1.28 

1.11 

1.06 


0.  73X10' 
.  73 
.74 
.77 
.88 
1.20 
2.  12 
5.  97 
21.0 


1.06X10' 
1.06 
1.08 
1.  11 
1.28 
1.74 
3.07 
8.66 
30.4 


0.515X10' 

.519 

.531 

.550 

.638 

.696 

.638 

.550 

.515 


0.747X10' 

.753 

.770 

.798 

.926 

1.01 

.926 

.798 

.747 


Coupling 

coefficients 

V 

cxs 

cys 

0 

0 

2.56 

.053 

3.00 

.  109 

2.62 

.  170 

1.44 

.410 

.793 

.793 

.410 

1.44 

2.62 

0 

L  Length  of  test  section  Iq  plus  tapered  portion  of  end  tabs 
i-j-  Length  of  tapered  portion  of  Micarta  end  tabs 

Section  just  beyond  tapered  portion  of  end  tabs  and  site  of  top  gages 

Figure  1.  -  Grid  for  finite-element  analysis  of  Mod  I/E  specimens.  (Top  gages  Ixated  at  nodes  74  and 
77;  midpoint  gages  located  at  nodes  326,  329,  and  332.  All  dimensions  shown  are  relative. ) 
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Figure  2.  -  Axial  stress  variation  at  tab  ends  for  Mod  I/E  specimens  for  several  load 
angles  (finite-element  analysis  using  the  experimental  fracture  load  for  each 
specimen). 


Load 


Figure  3.  -  Axial  strain  variation  at  tab  ends  for  Mod  I/E  specimens  for 
several  load  angles  (finite-element  analysis  using  the  experimental 
fracture  load  for  each  specimen). 
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Distance  across  specimen,  percent 


Figure  4.  -  Axial  stress  variation  at  midlengths  for  Mod  I/E  specimens  for  several 
load  angles  (finite-element  analysis  using  the  experimental  fracture  load  for 
each  specimen). 
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Figure  5.  -  Axial  strain  variation  at  midlengths  for  Mod  I/E  specimens 
for  several  load  angles  (finite-element  analysis  using  the  experi¬ 
mental  fracture  load  for  each  specimen). 
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deg 
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Figure  6.  -  Comparison  at  fracture  load  of  predicted  and  measured  axial 
strains  near  end  tabs  for  off-axis  specimens  for  Mod  I/E. 
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deg 

O  10^ 
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Distance  across  specimen,  percent 

Figure  7.  -  Comparison  at  fracture  load  of  predicted  and  measured  axial 
strains  at  midlength  of  off-axis  specimens  from  Mod  I/E. 


Figure  8.  -  NASTRAN  model  of  off-axis  specimen  (657  nodes;  576  CQUAD2  elements). 


(a)  Undeformed. 


(b)  Out-of-plane  bending  moment  (11.  30  N  -  m  (100  in  -  lb);  maximum  deflection,  4.  25  cm  (1. 672  in. )). 


(c)  Out-of-plane  twisting  moment  (11.  30  N  -  m  (100  in  -  lb);  maximum  deflection,  9. 29  cm  (3. 658  in. )). 

Figure  9.  -  NASTRAN  plots  of  the  10°  off-axis  specimen  showing  deformed  shapes  due  to  out-of-plane 
eccentricities  (Mod  I/E). 
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(a)  Undeformed. 


(b)  Out-of-plane  bending  moment  (11.  30  N*  m  (100  in*  lb);  maximum  deflection, 
16.35  cm  (6. 437  in. )). 


(c)  Out-of-plane  twisting  moment  (11. 30  N*  m  (100  in*  lb);  maximum  deflection, 

5.72  cm  (Z  250  in. )). 

Figure  10.  -  NASTRAN  plots  of  the  30®  off-axis  specimen  showing  deformed  shapes  due  to 
out-of-plane  eccentricities  (Mod  I/E). 
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Figure  12.  -  Out-of-plane  bending  and  twisting  effects  on  axial  strain  at 
midlength  of  10®  and  30®  off-axis  specimens  from  Mod  I/E  composite 
(11. 3-N-  m  (100-in  -  lb)  moments). 


Figure  13.  -  Schematic  depicting  Instrumentation  to  detect 
out -of -plane  eccentricities  during  testing  of  off-axis  fiber 
composites. 


Figure  14.  -  Comparison  of  finite-element  analysis  results  for  5®  off-axis 
specimen  (Mod  I/E)  showing  effects  of  specimen  thickness  variation. 


NASA-Lewis 


